Structures, stability, mechanical and electronic properties of a-boron and its twined 

brother a*-boron 
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The structures, stability, mechanical and electronic properties of a-boron and its twined brother 
a*-boron have been studied by first-principles calculations. Both a-boron and a*-boron consist 
of equivalent icosahedra B12 clusters in different connecting configurations of "3S-6D-3S" and 
"2S-6D-4S", respectively. The total energy calculations show that a*-boron is less stable than 
Q-boron but more favorable than /3-boron and 7-boron at zero pressure. Both a-boron and 
a*-boron are confirmed dynamically and mechanically stable. The mechanical and electronic 
properties of a-boron and a*-boron indicate that they are potential superhard semiconducting 
phases of element boron. 
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PACS numbers: 61.50.Ks, 61.66.-f, 71.20.Mq, 63.20. D- 

Boron is considered as the most mystic element due 
to its fascinating chemical complexity that leads to vast 
variety of polymorphs^"— . Among these polymorphs, 
rhombohedral a-boron and /3-boron (with 12 and 106 
atoms in unit cells, respectively) are well-known and 
almost degenerate in energy at ambient condition, 
a-boron is more favorable than /3-boron at high pressure 
and /3-boron is more stable in high temperature con- 
dition. In 2008, Oganov et al^^— discovered an ionic 
phase of elemental boron (orthorhombic 7-boron with 28 
atoms in its unit cell), which consists of icosahedral B12 
clusters and interstitial B2 pairs acting as anions and 
cations, respectively, in a NaCl-type arrangement. This 
new phase of boron is stable between 19 and 89 GPa but 
quenchable at ambient conditions, which confirms the 
previous discovery of a new crystal of boroni and thus 
provides us the missing piece of a puzzle of the phase 
diagram of boron-. 

Born atoms in its element crystals prefer forming 
icosahedral B12 cluster, where each boron atom connects 
five neighbors with intra-Bi2 B-B bonds. In a-boron, 
equivalent icosahedral B12 clusters connect to each 
other directly with inter-Bi2 B-B bonds and each of 
them has 12 equivalent neighbouring B12 clusters. In 
7-boron, each icosahedral B12 connects to ten equivalent 
neighbouring B12 clusters directly with inter-Bi2 B-B 
bonds and two secondly adjacent B12 clusters indirectly 
by interstitial B2 pairs. The connectivity between icosa- 
hedral B12 clusters and interstitial atoms in /3-boron 
becomes more complex. 

a-boron is not the only phase of boron consisting of 
pure icosahedral B12 clusters (without any interstitial 
atoms like those in /3-boron and 7-boron) . Very recently, 
Pickard et al. obtainedii a new meta-stable phase 
of boron by an ab initio random structure searching 
method. This new phase belongs to Cmcm space group 



and can be considered as a twined polymorph of a-boron 
with new connectivity (different from that in a-boron 
boron) between icosahedral B12 clusters which is 0.01 
eV/atom less stable than a-boron and more favorable 
than /3-boron. With remarkable energetic stability, such 
a phase of boron is a potential superhard material. 
But there is no further discussion about the structure, 
dynamical stability, mechanical and electronic properties 
of this new phase. In this paper, based on the density 
functional theory (DFT) method, we systematically 
study the structure, dynamical and mechanical stability, 
mechanical and electronic properties of this new boron 
phase (we call it a*-boron in the following) and compare 
them with those of a-boron. 



I. COMPUTATIONAL DETAILS 

All calculations are carried out using the density func- 
tional theory (DFT) with general gradient approximation 
(GGA)-- as implemented in Vienna ab initio simulation 
package (VASP)i^'i^. The interactions between nucleus 
and the valence electrons are described by the projector 
augmented wave (PAW) methodi^ii^. A plane-wave ba- 
sis with a cutoff energy of 500 eV is used to expand the 
wave functions of all the systems considered in present 
work. The Brillouin Zone (BZ) sample meshes for all 
systems are set to be denser enough (less that 0.2 1/A) 
in our calculations. Crystal lattices and atoms positions 
for every systems are fully optimized (under different ex- 
ternal pressure) up to the residual force on each atom 
less than 0.005 eV/A through the conjugate-gradient al- 
gorithm. The vibrational properties of a-boron and a*- 
boron are studied by using the phonon package^^ with the 
forces calculated from VASP. For the mechanical proper- 
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FIG. 1: Primitive view and crystal view of a-boron (a) ; Sketch 
maps showing the structural characters of a-boron (b) and 
Q:*-boron (c); Primitive view and crystal view of a*-boron 
(d); Simulated X-ray (A= 0.31851 A) diffraction patterns of 
a-boron and a*-boron (e). In figures the blue and brown balls 
indicate the S-type and D-type boron atoms, respectively. 

ties calculations, different numbers of deformed a-boron, 
a*-boron and 7-boron are considered according to their 
crystal symmetries. 



II. RESULTS AND DISCUSSIONS 

A. Structures 

The optimized structures of a-boron and a*-boron in 
both primitive cell (top) and crystal cell (bottom) are 
shown in Fig. [T] (a) and (d), respectively. The hexago- 
nal crystal cell of a-boron belongs to R-3m (166) space 
group. Under zero pressure, its lattice constants are 
a=b=4.9 A and c=12.55 A. Two inequivalent boron 
atoms Bl and B2 in a-boron locate at positions of (0.237, 
0.119, 0.108) and (0.197, 0.803, 0.024), respectively. In a- 
boron, each Bl atom possesses five intra-Bi2 B-B bonds 
and single (S-) inter-Bi2 B-B bond (S-type, colored in 
blue). Each B2 atom possesses five intra-Bi2 B-B bonds 
and double (D-) inter-Bi2 B-B bonds (D-type, colored in 
brown). The orthorhombic a*-boron belongs to Cmcm 
(63) space group. Its equilibrium lattice parameters are 
a=4.88 A, b=8.85 A and c= 8.06 A, respectively. There 
are five inequivalent boron atoms Bl, B2, B3, B4 and 
B5 in a*-boron. They locate at positions of (0.667, 
0.334, 0.75), (0.818, 0.507, 0.75), (1.000, 0.563, 0.924), 
(0.797, 0.169, 0.639) and (0.500, 0.265, 0.568), respec- 
tively. Similar to those in a-boron, the five inequivalent 



boron atoms in a*-boron can be divided into two types: 
S-type where the three boron atoms (Bl, B2 and B3) 
colored in blue have five intra-Bi2 B-B bonds and sin- 
gle (S-) inter-Bi2 B-B bonds and D-type where the two 
brown boron atoms (B4 and B5) possess five intra-Bi2 
B-B bonds and double (D-) inter-Bi2 B-B bonds. Al- 
though a-boron and a*-boron contain different numbers 
of inequivalent boron atoms, both structures contain only 
one inequivalent icosahedral B12 cluster. That is the rea- 
son we refer them as twined structures of element boron. 

In both a-boron and a*-boron, each icosahedral B12 
cluster possesses twelve neighboring B12 clusters dis- 
tributing in different stacking manners. In a-boron, 
twelve neighboring B12 clusters distribute around the 
center B12 cluster with a "3S-6D-3S" connecting config- 
uration. As indicated in Fig [1] (b), three top-position 
B12 clusters (indicated as 7, 8 and 9) connect to the cen- 
ter B12 cluster with single inter-Bi2 B-B bonds formed 
by S-type boron atoms, six middle-position B12 clusters 
(indicated as 1, 2, 3, 4, 5 and 6) connect to the cen- 
ter B12 clusters with double inter-Bi2 B-B bonds formed 
by D-type boron atoms and three bottom-position B12 
clusters (indicated as 10 ,11 and 12) connect to the cen- 
ter B12 cluster with single inter-Bi2 B-B bonds formed 
by S-type boron atoms, forming a "3S-6D-3S" configu- 
ration. In a*-boron, the twelve neighboring B12 clusters 
distribute around the center B12 cluster with a "2S-6D- 
4S" connecting configuration as indicated in Fig [1] (c). 
Two top-position B12 clusters (indicated as 7 and 8) con- 
nect to the center B12 cluster with single inter-Bi2 B-B 
bonds formed by S-type boron atoms, six middle-position 
B12 clusters (indicated as 1, 2, 3, 4, 5 and 6) connect to 
the center B12 cluster with double inter-Bi2 B-B bonds 
formed by D-type boron atoms and four bottom-position 
B12 clusters (9, 10, 11 and 12) connect to the center B12 
cluster with inter-Bi2 B-B bonds formed by S-type boron 
atoms. In both a-boron and a*-boron, S-type atoms con- 
nect to D-type atoms only with intra-Bi2 B-B bonds and 
never meet each other in different B12 clusters. 

Both a-boron and a*-boron contain equivalent icosa- 
hedral B12 clusters in their crystal structure. The dif- 
ferent connecting configurations of "3S-6D-3S" and "2S- 
6D-4S" for a-boron and a*-boron result in different space 
group of R-3m (166) and Cmcm (63), respectively. The 
different structures of a-boron and a*-boron also result 
in different X-ray diffraction patterns. Fig. [T] (e) shows 
the X-ray diffraction patterns of a-boron and a*-boron. 
We can see that almost all the diffraction peaks of a- 
boron can be found in the XRD of a*-boron with re- 
duced intensities, such as the peaks at 4.5°, 7.2° and 8.6°. 
In addition, a*-boron possesses characteristic diffraction 
peaks forbidden in a-boron, which locate at 4.7°, 6.2°, 
8.0°, 9.9°, 12.09°, 19.05° and 23.33°. These character- 
istic peaks are helpful for experimental identification of 
a*-boron. 
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TABLE I: The calculated values of elastic constants dj (GPa), bulk modulus B (GPa), shear modulus G (GPa), Vicker's 
hardness H„ (GPa) of the Q-boron, a*-boron and 7-boron. 
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FIG. 2: The enthalpy curves for Qf*-boron, 7-boron, /3-boron 
and a-Ga-boron as functions of pressure relative to that of a- 
boron (a); The simulated phonon band structures of a-boron 
(b) and Q*-boron (c). 



B. stability 

In this part, we discuss the energetic and dynamical 
stability of a-boron and a*-boron. In Fig. [2] (a) we 
show the enthalpy curves of a^-boron, 7-boron, /3-boron 
and a-Ga-boron as functions of pressures relative to 
that of a-boron. We can see that a*-boron possesses 
remarkable energetic stability comparable to that of 
its twined brother a-boron at any pressure range. Its 
enthalpy is only about 0.01 eV/atom above that of a- 
boron. At zero pressure, both a-boron and a*-boron are 
more favorable than 7-boron, /3-boron and a-Ga-boron. 
When the external pressure is increased, they transform 
to 7-boron at 20 GPa and then to a-Ga-boron at 94 
GPa. Our results are in good agreement with previous 
first- principles calculations^. 

The remarkable stability indicates that a*-boron has 
high probability to coexist with a-boron as the ground 
state of element boron in the pressure range from 0-21 
GPa. But how about the possibility? Is it a dynamically 
stable phase? We need to know the vibrational prop- 
erties of such a new crystalline structure of boron to 
answer the question. We investigate the phonon prop- 
erties of a-boron and a*-boron by phonon package^^ 
with interatomic forces calculated from vaspi^ii^. The 



calculated phonon band structures of a-boron and 
a*-boron are shown in Fig. [5] (b) and (c), respectively. 
We can see that there are no imaginary frequencies 
in both a-boron and a*-boron, indicating that they 
are dynamically stable. Above discussions about the 
energetic and dynamical stability of a*-boron indicate 
that it is a promising meta-stable phase of element boron. 



C. Mechanical properties 

The elastic constants of a-boron and a*-boron are 
calculated as the second-order coefficient in the poly- 
nomial function of distortion parameter S used to fit 
their total energies according to the Hooke's law. In 
view of their differences in crystal symmetries, six and 
nine groups of deformations are applied on a-boron 
and a*-boron (7-boron) in our calculations, respec- 
tively. The elastic constants of a-boron and a*-boron 
are summarized in Table I with those of 7-boron used 
for comparison. Our calculated elastic constants for a- 
boron and 7-boron are in good agreement with those 
in previous first-principle calculations^!^. Based on 
these results, we can see that the six independent elas- 
tic constants of a-boron satisfy the corresponding me- 
chanical stability criteria^®"^° for a hexagonal crystal 
with C66=(Cii-Ci2)/2>0, Cii+Ci2-2(Ci3)VC33>0, and 
C44>0. The orthorhombic a*-boron possesses nine 
independent elastic constants. As listed in Table I, 
these elastic constants also satisfy the corresponding me- 
chanical stability criteriai^"— for a orthorhombic crys- 
tal with Cii-|-Ci2-C2l>0, Cii-hC33-2Ci3>0, C22 + C33- 
2C23>0, Cii>0, C22>0, C33>0, C44>0, C55>0, C66>0, 

and Cii-HC22+C33+2Ci2+2Ci3+2C23>0. 

The bulk modulus (B) and shear modulus (G) are 
evaluated according to Hill's formula^^ based on the cal- 
culated elastic constants. The calculated B and G for 
a*-boron are 205.58 GPa and 190.91 GPa, respectively, 
which are very close to those of a-boron. To further an- 
alyze the hardness of a-boron and a*-boron, we adopt 
the recently introduced empirical scheme^^ to evaluate 
their Vicker's hardness (H^) determined their B and G as 
H.„=2(GVB2)°-^^'5-3. The calculated values of Vicker's 
hardness for a-boron and a*-boron are 38.96 GPa and 
36.60 GPa, respectively. These values indicate that both 
a-boron and a*-boron are potential superhard phases of 
element boron. 



4 




S R 



PDOS 



PDOS 



FIG. 3: The band structures of a-boron (a) and corresponding 
local density of state projected on the S-type (b) and the D- 
type (c) boron atoms. The band structures of a*-boron (d) 
and corresponding local density of state projected on the S- 
type (e) and the D-type (f) boron atoms. 
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boron, respectively. From the band structures, we can 
see that both of them are indirect-band-gap semiconduc- 
tors. Differently, the valence band maximum (VBM) and 
conduction band minimum (CBM) in a-boron locate at 
T-point (0.5, 0.5, 0.5) and T-point (0.0, 0.0, 0.0), re- 
spectively, separated by a forbidden band with a gap of 
1.497 eV. But the VBM and CBM in a*-boron locate at 
r-point (0.0, 0.0, 0.0) and T-point (0.5, 0.5, 0.5), respec- 
tively, forming a indirect-band gap of 2.004 eV. Based 
on the local density of states (LDOS) projected on S- 
type and D-type boron atoms as shown in Fig. [31 we can 
see that the VBM states in both a-boron and a*-boron 
are contributed from the D-type boron atoms and their 
CBM states are derived from S-type boron atoms. In 
both a-boron and a*-boron the whole valence band is 
mainly contributed by the S-type boron atoms and the 
whole conduction band is mainly contributed by the D- 
type boron atoms. These results indicate that the valence 
electrons density mainly comes from the S-type boron 
atoms and the hole density mainly comes from the D- 
type atoms. 

The corresponding charge difference density (CDD) 
of a-boron and a*-boron is shown in Fig. 2] to confirm 
and understand the electronic properties discussed above. 
From the CDD of a-boron shown in the left in Fig. HJ we 
can see that the D-type boron atoms (brown balls) loss 
electrons and the S-type boron atoms (blue balls) gain 
electrons. The CDD of a*-boron is shown in the right 
panel in Fig. Similar to that of a-boron, the D-type 
boron atoms (brown balls) in a*-boron loss electrons and 
the S-type boron atoms (blue balls) gain electrons. These 
CDD results reveal that the bonding electrons in both a- 
boron and a*-boron prefer to distribute on the inter-Bi2 
B-B bonds formed by the S-type boron atoms than those 
formed by the D-type boron atoms and the intra-Bi2 B-B 
bonds. 



III. CONCLUSION 



FIG. 4: The charge difference density for Q-boron (left) and 
a*-boron (right). Electrons transfer from the green- blue area 
to the red-yellow area. 



D. Electronic properties 

Both a-boron and a*-boron contain pure icosahe- 
dral Bi2 clusters with different connecting configurations. 
Such a structural difference may result in different elec- 
tronic properties in a-boron and a*-boron. In Fig. |3] (a) 
and (d), we show the band structures of a-boron and a*- 



We have studied the structure, stability, mechanical 
and electronic properties of a new orthorhombic boron 
phase a*-boron and compared them with those of its 
twined brother a-boron. Both a-boron and a*-boron 
consist of equivalent icosahedra B12 clusters with 
different connecting configurations of "3S-6D-3S" and 
"2S-6D-4S", respectively. Our results show that a*- 
boron is a dynamically and mechanically viable phase of 
element boron. The remarkable energetic stability and 
considerable hardness of 36.8 GPa of a*-boron indicate 
that it is a promising superhard material with potential 
applications in industry. 
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